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Abstract 

In this paper, we study the response of 2D framed structures made of rectangular cells, to the sudden 
removal of columns. We employ a simulation algorithm based on the Discrete Element Method, where 
the structural elements are represented by elasto-plastic Euler Bernoulli beams with elongation-rotation 
failure threshold. The effect of structural cell slenderness and of topological hierarchy on the dynamic 
residual strength after damage Ri is investigated. Topologically hierarchical frames have a primary 
structure made of few massive elements, while homogeneous frames are made of many thin elements. We 
also show how Ri depends on the activated collapse mechanisms, which are determined by the mechanical 
hierarchy between beams and columns, i.e. by their relative strength and stiffness. Finally, principles of 
robustness-oriented capacity design which seem to be in contrast to the conventional anti-seismic capacity 
design are addressed. 
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Introduction 


Since many decades, design codes ensure a very low probability that a building collapses under ordinary 
loads, like self weight, dead and live service load, or snow. Nevertheless buildings still do collapse, from time 
to time. An extremely small fraction of collapses originates from unlikely combinations of intense ordinary 
load with very poor strength of the building. The majority of structural collapses are due to accidental events 
that are not considered in standard design. Examples of such events are: gross design or construction errors, 
irresponsible disregard of rules or design prescriptions, and several rare load scenarios like e.g. earthquakes, 
fire, floods, settlements, impacts, or explosions [AlexandCT 2004 . Accidental events have low probability 
of occurrence, but high potential negative consequences. Since risk is a combination of probability and 


consequences, the risk related to accidental events is generally significant BS Eurocode 1, 2004 


In 1968 a gas explosion provoked the partial collapse of the Ronan Point building in London. This event 
highlighted for the first time the urgency for robust structures, enduring safety in extraordinary scenarios 

Since then, interes was driven by striking catastrophic collapses Val and Val[ 


Pearson and Delatte, 2005 


2006|, until in 2001 the tragic collapse of the World Trade Center renewed the attention to the topic (see 
e.g. [Bazant and Zhou 2002 and Cherepanov and Esparragoza 2007| ). The last decades, several design 
rules aimed at improving structural robustness have been developed (see e.g. Masoero 2010|) . 

Accidental events can be classified into identified and unidentified IBS' Eurocode 1 


2004 DoD 2005 


Identified events are statistically characterizable in terms of intensity and frequency of occurrence. Examples 
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are earthquakes, fire not fueled by external sources, gas explosions, and unintentional impacts by ordinary 
vehicles, airplanes, trains, or boats. Specific design rules and even entire codes are devoted to specific 
identified accidental events. Unidentified events comprise a wide variety of incidents whose intensity and 
frequency of occurrence can not be described statistically, e.g. terrorist attacks or gross errors. 

The risk related to unidentified accidental events can be mitigated both by structural and nonstructural 
measures Gulvanessian and Vrouwenvelder 2006 . Nonstructural measures such as barriers and monitoring 


can reduce the probability that an accidental event affects the structural integrity, others like a wise dis¬ 
tribution of plants and facilities can minimize the negative consequences of eventual collapses. Otherwise, 
structural measures can improve local resistance of structural elements to direct damage, e.g. the design of 
key elements for intense local load BS Eurocode 1, 2004 , or the application of the Enhanced Local Resis¬ 
tance method DoD, 2005 . Structural measures can also provide progressive collapse resistance, i.e. prevent 


spreading of local direct damage inside the structure to an extent that is disproportioned with respect to 
the initial event. Usual strategies to improve progressive collapse resistance are compartmentalization of 

and delocalization of stress after local damage. Stress delocalization can be 


structures 


Starossek 2006 


obtained exploiting redundancy, plastic stress redistributions (Masoero, Wittel et ah, 2010), ties [Alexander 


2004 , and moment resisting connections Hamburger and Whittaker, 2004 Vlassis et al., 2008] . 


Nowadays several design codes employ the conventinal Alternate Load Path Method (ALPM) to evaluate 
progressive collapse resistance, e.g. GSA 2003 and DoD 2005 . The method consists in removing one key 


element, generally a column or a wall, and measuring the extent of subsequent collapse. If the final collapse 
is unacceptably wide, some of the previously listed measures have to be employed. Hence structures are first 
designed and subsequently tested to be robust - they are not conceived a priori. This course of action excludes 
optimizations of the basic structural topology and geometry, that actually play a key role in the response to 
local damage, considering as an example the very different behavior of redundant and statically determined 
structures. Anti-seismic design BS Eurocode 8 2004 already contains some prescriptions that should 


be considered before starting a new design, e.g. geometric regularity on the horizontal and on the vertical 
planes. Furthermore, anti-seismic capacity design requires a hierarchy of the structural elements ensuring that 
earthquakes can only provoke ductile collapse of the horizontal beams, while failure of columns and brittle 
ruptures due to shear are inhibited. Differently, for what concerns progressive collapse resistance, optimal 
overall geometric features are not known, except for the concepts of redundancy and compartmentalization. 
Furthermore, the idea of hierarchically maximizing progressive collapse resistance is completely absent. 

In this paper, we make a first step to cover this deficiency, showing that progressive collapse resistance 
can be improved by hierarchy in the overall geometry {topological hierarchy) and in the relative strength and 
stiffness of horizontal and vertical structural elements {mechanical hierarchy). Our approach incorporates 
the simulation of progressive collapse of regular 2D frames made of reinforced concrete (RG) subjected to 
the sudden removal of structural elements, following the ALPM framework. We first describe the analyzed 
frame structures and briefly sketch the approach that is based on the Discrete Element Method (DEM). 
After the model description, we present the results of the simulations, with focus on the effect of geometry 
and hierarchy on the activated collapse mechanisms and, consequently, on progressive collapse resistance. 


Hierarchical structures and damage 

We consider two representative sets of regular 2D framed structures in Fig. [^a. Each set consists of three 
frames with identical total width Ltot and different topological hierarchical level 1/n, where is the number 
of structural cells in a frame. The horizontal beams, excluded those of the secondary structure, carry a 
uniform load per unit length qext- The frames are made of RG with typical mechanical parameters of 
concrete and steel, as shown in Table [l] The total height Htot of the structure is kept constant, and two 
different height-bay aspect ratios A = H/L of the structural cells are considered (A = 0.75 and A = 1.33). 

There exist several ways of introducing hierarchy into the topology of framed structures; here we call a 
structure “hierarchical” if it has a primary structure, made of few massive structural elements, that supports 
a secondary one. The latter defines the living space and has negligible stiffness and strength compared to the 
primary structure. The frames with n = 2 and n = 5 can be seen as reorganizations of those with n = 11. 
In detail, each column of the frames with n = 5 corresponds to two columns of the frames with n = 11, and 
the same is valid for the beams, disregarding the first floor beam of the frames with n = 11, which is simply 
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Figure 1: a) Map of the studied frames. The cohesion of the elements inside the dotted damage area is 
suddenly removed to represent the initial accidental damage, b) Cross sections and rebars of beams (above) 
and columns (below). 


deleted (see Fig. [^a). Analogously, the geometry of the frames with n = 2 can be obtained starting from 
the frames with n = 5. 

The cross sections of columns are square (see Fig. Sb), with edges he = be proportional to H with 
factor Ac = 1/10. The beams have rectangular cross section whose height ht is proportional to L with factor 
Xb = 1/10, and whose base bt is proportional to hb with aspect ratio 6b = 2/3. The reinforcement is arranged 
as shown in Fig. 0b, with area Ag proportional to the area of the cross section by factor ps^e = 0.0226 for 
the columns (i.e. 8(/18 when n=ll), and ps^b = 0.0029 for the beams (i.e. 4(/)14 when n=ll). 

The damage areas, (dotted in Fig. [^a) contain the structural elements that are suddenly removed to 
represent an accidental damage event, following the ALPM framework. The damage is identical for frames 
with same A, and is defined by the breakdown of one third of the columns on a horizontal line. The columns 
and beams removed from frames with n = 11 correspond to the structural elements removed from frames 
with n = 5 and n = 2. This kind of damage is employed to represent accidental events with a given 


Parameter 

Symbol 

U nits 

Value 

Reinforced concrete 

Specific weight 

IRC 

kg/m^ 

2500 

Young modulus 

Ee 

N/m2 

30-10® 

Compressive strength (high) 

fc 

N/m2 

35-10® 

Compressive low (low) 

fc 

N/m2 

0.35-10® 

Ultimate shortening 


- 

0.0035 

Steel 

Young’s modulus 

Eg 

N/m^ 

200-10® 

Yield stress 

fy 

N/m2 

440-10® 

Ultimate strain 

^U,S 

- 

0.05 


Table 1: Mechanical properties of reinforced concrete and steel. 
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amount of destructive energy or spatial extent, like explosions or impacts. In this work we do not explicitly 
simulate very local damage events like gross errors, which would be better represented by the removal of 
single elements. Nevertheless, we will generalize our results to consider also localized damage events. 


DEM model 


We employ the Discrete Element Method (DEM) to simulate the dynamics after a sudden damage Poschel 
and Schwager 2005 Carmona et al. 2008] . DEM is based on a Lagrangian framework, where the structure 
is meshed by massive elements interacting through force potentials. The equations of motion are directly 
integrated, in our case using a 5**' order Gear predictor-corrector scheme, with time increments between 
10“®s and 10“®s (see Masoero, Wittel et ah, 2010). DEM is an equivalent formulation to Finite Elements, 
converging to the same numerical solution of the dynamics if identical force-displacement laws are imple¬ 
mented. A detailed description of the algorithm for 3D systems can be found in (Masoero, Wittel et ah, 
2010); [M asoero[ |2010| , together with a discussion on the applicability. For this work, the code was restricted 
to 2D by allowing only two displacements and one rotation in the vertical plane. In (Masoero, Vallini et ah, 
2010), the DEM model is tested against dynamic energy-based collapse analyses of a continuous horizontal 
beam suddenly losing a support. In the appendix, we compare our DEM results to experimental observa¬ 
tions of a 2D frame undergoing quasi-static column removal Yi et al. 2008 . To the best of our knowlede. 


literature still lacks on experiments of dynamic collapse of framed structures due to accidental damage. 

In the following we will review only the essentials of our model, focusing on the details tha are relevant 
for the application to 2D frames. We assume simplified force-displacement laws for the beam element and 
for the Hertzian contacts. Predicting collapse of real structures would require more specialized interaction as 
compared to here, for example using the fiber approach for the cross sections. By contrast, we are interested 
in fundamental mechanisms of damage propagation within complex structural systems. In this research 
perspective, and according to a basic principle of Statistical Mechanics, minimizing the complexity of local 
interactions improves the interpretation of the systemic response. Despite the strong assumptions, in the 
Appendix we show that our model can match reasonably well with with experimental observations. 


Structural representation 

In a first step, the structure needs to be assembled by discrete elements and beams. Fig.j^a shows the four 
types of Spherical Discrete Elements (SDE) that we employed. Columns and beams are made of 9 SDEs, 
respectively with diameter = 0.8/ic and Df, = 0.8/ib, slightly smaller than the distance between them to 
prevent contact form occurring before local rupture. Constrained SDEs and connection SDEs have same 
diameter Dc as column SDEs. The constrained SDEs are clamped to a plane that represents the ground by 
means of the Hertzian contact model, discussed further in this section. 




1 Connection: 

! Column! : 

: SDEs : 

; SDEs i . 



' Beam 

\ Constraint ' 

i SDEs 

: SDEs 




Ground plane 

(a) 



Figure 2: a) Typical DEM mesh of a 2D frame. The centers of the spheres are connected by b) Euler-Bernoulli 
beam elements, shown in undeformed (above) and deformed state (below). 
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Euler-Bernoulli beam elements 


Pairs of SDEs are connected by Euler-Bernoulli beam elements - (EBE) that, when deformed, transmit 
forces and moments to their edge nodes, locally labeled 0 and 1 (see Fig. §b). The mass M of an SDE is 
defined on the basis of the EBEs connected to it. Namely M = X) 1/2 where e labels the generic 

e 

EBE connected to SDE, and ^4® is the cross sectional area of the structural element corresponding to the 
EBE. The external load q^xt is introduced adding a mass qextL^lg to the beam SDEs. q^xt is not treated 
directly as a force to avoid downward accelerations of the SDEs greater than gravity g during free fall. 

For sufficiently small deformations, the EBEs are linear elastic and exert a force N proportional to 
the elongation e and directed along the 01 segment, a shear force T proportional to the sum of the nodal 
rotations 6 ((^o + </^i)) and a bending moment B proportional to the nodal effective rotations, defined as 
= 4(/Jo + 2ipi, and = 2ipQ + Furthermore, we introduce damping by forces and moments 

directed opposite to N, T, and B, and proportional to the time derivative of e with factor 7 ^ = lOONs/m, 
and of (po,ipi with factor jb = lONms. Geometric nonlinearity due to large displacements is considered by 
referring rotations and elongation to the 01 segment. In the small deformations regime of our simulations, 
N is with good approximation equal to the axial force inside the EBE, and thus perpendicular to T. 

If N overcomes a yield threshold in tension Ny or under compression N^y, the ideally plastic regime is 
entered and plastic axial strain is applied to maintain N = Ny or N = N^y. Neglecting the contributions 
of concrete in tension and of steel in compression, we set the yield thresholds in terms of N and e to: 


Ny - Asfy 

N,y = (A" 


r - and 

^ ~ A-Er " " ’ 


Er 


As)fc 


N, 


cy 


{A--As)E, 


A 

Ec ■ 


( 1 ) 

( 2 ) 


Ideally plastic regime in bending is entered when \B\ > By. We obtain the bending yield threshold By 
and the corresponding yielding effective rotation q^y^^, neglecting the strength contribution of concrete and 
assuming a lever arm between upper and lower reinforcement equal to the height h of the cross section: 


By = tsPsA^fyh + ABy und 




Ey 

EJ^ 


(3) 


N is the cross sectional moment of inertia of the EBE, and G is the fraction of reinforcement in tension (3/8 
for columns and 1/2 for beams, as in Fig. Sb). ABy considers the beneficial compression effect compression 
in the EBE. We set ABy assuming bending carried by the reinforcement alone, and that the strain Eg {ABy) 
in the reinforcement put under tension by ABy equals the compressive strain Eg {N) due to <0, namely: 


ffs (iV) = £« (ABy) 


N _ ABy 

~ tgpgA<^hEg 


(4) 


In this way, eventual tension > 0 inside the EBE produces negative ABy, and thus reduces By. When 
yielding in bending occurs, plastic rotations are added at the edge nodes of the EBE. If only \Bi\, with 
j = 0,1, is greater than By, then only ip^ is applied to restore \Bi\ = By. Differently, if both |i3o| and |i?i| 
are greater than By, both ip^ ^nd are applied to restore |i3o| = |i?i| = By. For the sake of simplicity, 
we assume yielding in bending uncoupled from yielding in axial direction. Furthermore, we neglect yielding 
due to shear because small plastic deformations are generally associated with shear. 

We consider an EBE failed when excessive e^ and pE are cumulated. For this purpose, the coupled 
breaking criterion: 


rPl 


max 


pi 

Ei 


i^th ^y) ^ \ Pth 


-pi 




max 


P‘, 


pi 


Pth 


>1 if > 0 and 

>1 if £P' < 0 , 


(5) 

( 6 ) 
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is employed. Sth, Scth, and ipth are the maximum allowed plastic elongation, shortening, and rotation in 
uncoupled conditions. We consider high plastic capacity of the structural elements setting eth = 

£cth = ‘2eu,c, and tpth = 0.2rad (see Table [^. Failed EBEs are instantly removed from the system. We 
neglect ruptures due to shear assuming that, in agreement with a basic principle of capacity design, a 
sufficient amount of bracings ensures the necessary shear strength. 


Inter-sphere contact 


The Hertzian contact model is employed for the SDEs to consider collisions between structural elements. 
The model consist of repulsive forces between partially overlapping SDEs, damped by additional forces 
proportional and opposite to the overlapping velocity. We also set tangential forces that simulate static and 
dynamic friction, as well as damping moments opposed to the relative rolling velocity. A similar Hertzian 
contact model is also employed for SDEs colliding with the ground plane. In the following simulations 
we employ contact parameters that can be found in Masoero 2010 . We do not transcribe them because 


impacts do not affect significantly the collapse mechanisms sudied here. Nevertheless, in general simulation 
algorithms for progressive collapse should consider impacts, because initial damage located at upper stories 
generates falling debris, and because impacts can drive the transition from partial to total collapse (see 
(Masoero, Wittel et ah, 2010) and Bazant and Zhou 2002|). In granular dynamics, the contact parameters 


are generally set referring to the material of the grains [Poschel and Schwager 2005 . In our model the 


SDEs represent large heterogeneous portions of structural elements, for which there are not conventionally 
defined contact parameters so far. We emply parameters yielding a qualitatively realistic dynamics (e.g. the 
elements do not rebound or pass through each other), and chosen from sets of possible one that were defined 
through preliminary studies. Such studies also indicated that the collapse loads of a beam due to debris 
impact varies of less than 15% upon orders of magnitude change in the contact parameters. 


Simulating progressive collapse 


The simulations are organized into two steps: first the structure is equilibrated under the effect of q^xt ^^nd 
gravity, then the EBEs inside the damage area are suddenly removed, and the subsequent dynamic response 
is simulated. Our aim is to quantify three collapse loads: 

• q^: maximum static load that the intact structure can carry; 

• qc'. minimum critical load that causes dynamic collapse after damage. Applied statically to the intact 
structure first, it is then kept constant during the post-damage dynamic response. 

• qp^t'- minimum load corresponding to total collapse after damage. By definition, q^ < qp^t < Qu¬ 


in our DEM model we do not have a straightforward unique measure of load, because the mass of the SDEs 
depends on the external load q^xt and on the self weight of the structural elements. The mass of the beam 
SDEs effectively acts as a distributed horizontal load. On the other hand, the columns at each story transmit 
vertical concentrated forces either to other columns at a lower story, or to the horizontal transfer beam over 
the damage area. Therefore we introduce a load measure that we call equivalent load q^q, applied to the 
massless structure and analytically related to the geometry, the mass, and the activated collapse mechanism 
of the frames in the simulations. Namely, qgq is defined to produce the same static effect as the various 
masses and concentrated forces of the simulation frames, at the critical points where collapse is triggered. 


The derivation of the analytical expressions used in this work is shown in Masoero 2010 


For each analyzed structure, we first apply the entire structural mass. In a subsequent step, the external 
load qext is increased until the intact structure collapses in static conditions. The collapse mechanism 
indicates what equivalent load expression should be used to compute Then we slightly decrease qext, 
equilibrate, introduce the damage, and calculate whether dynamic progressive collapse is triggered and to 
what an extent. Performing several simulations with progressively smaller qext, the final extent of collapse 
changes from total to partial, and we employ again an adequate equivalent load to compute qp^f If the 
structure collapses even when qext is reduced to zero, we start reducing the specific weight of the structural 
elements, i.e. the structural mass. When dynamic collapse does not occur anymore, an adequate equivalent 
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load provides qc- Once we obtain the collapse loads, we estimate the progressive collapse resistance referring 
to the residual strength fraction Ri = Qc/qi- Actually, progressive collapse resistance is more directly related 
to qc, but the advantage of i?i is that it can not be improved by simply strengthening the structural elements, 
which would increase both qc and Robustness-oriented structural optimization is required to increase 
i?i, which therefore is a good indicator to compare different structural solutions. 


Bending collapse 

In our model, the bending yield threshold By does not depend on the strength of concrete fc- Therefore, 
setting the high value fc = 35N/mm^, the mainly compressed columns get much stronger than the horizontal 
beams, that fail in bending (see Figs. 3j4|. The resulting collapse mechanisms resemble triple-hinge and four 
hinges mechanisms, reflecting the large plastic capacity of the structural elements. 
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Figure 3: Static bending collapse mechanism before damage. Time 
of an EBE. 
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Figure 4: Dynamic bending collapse mechanism after damage. Time / = Os corresponds to the application 
of the initial damage. 


If the initial damage triggers a bending mechanism, frames with n = 2 undergo total collapse, while 
frames with lower hierarchical level 1/n initially suffer only a local collapse (see Fig. [^. The local collapse 
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can nevertheless evolve to total collapse, if high applied load and plastic capacity cause the falling central 
part of the structure to dynamically drag down the lateral portions (Masoero, Wittel et ah, 2010). 




Figure 5: Snapshots of partial ((7e9=13kN/m) and total (geq=26kN/m) bending collapse after damage of a 
frame with very strong columns, n=ll, and A = 1.33. Time t = Os corresponds to the application of the 
initial damage. 


The collapse loads, expressed in terms of equivalent loads qeq, are summarized in Fig. ^ as a function 
of the hierarchical level 1/n, for different slenderness of the structural cells A. In Fig. [6 superscript B 
indicates bending collapse mechanism. We employ equivalent loads referring to perfectly brittle or perfectly 
plastic bending failure (see the Appendix). The collapse loads decrease with A, i.e. a slender structure seems 
weaker, and increase with 1/n, i.e. hierarchical frames are stronger. The residual strength fraction Ri does 
not depend on A, while hierarchical structures with low n are more robust than homogeneous ones (see 
Fig.|§. In fact, the concentration of bending moment at the connection between a beam hanging above the 
damage area and the first intact column depends on the number of removed columns. In the simulations, 
we remove a constant fraction of one third of the columns on a horizontal line (see Fig. [^. Therefore 
homogeneous structures lose more columns and are less robust toward the bending collapse mechanisms. On 
the other hand, since the number of removed columns is decisive, we expect that the hierarchical level does 
not influence Ri toward bending collapse in case of single column removal. Finally we consider the 2D frame 
as part of a regular 3D structure and divide the collapse loads in Fig. by L, i.e.by the tributary length of 
the beams in the direction perpendicular to the frame. In this way, collapse loads per unit area are obtained 
in Fig. showing that: A does not influence qc/L and q^/L] structures with slender cells are less likely to 
collapse entirely; q^/L is independent from the hierarchical level; qdL is proportional to 1/n. 

Pancake collapse 

Progressive compressive failure of the columns, also called pancake collapse, occurs when we set the com¬ 
pressive strength of concrete to a small value /c = 0.35N/mm^ (see Fig. |^. This choice is unphysical but 
allows us to separate the effect of strength reduction from that of stiffness reduction in the columns. More 
realistic scenarios would involve columns with small cross section and highly reinforced, tall beams. 

The columns immediately next to the damage area are the first to fail under compression, and then 
progressive collapse spreads horizontally to the outside. We employ equivalent loads qeq referring to the 
two limit cases of local and of global pancake collapse. Local pancake collapse occurs when the bending 
stiffness of the beams is very low and when the compressive failure of the columns is very brittle. In this 
case, the overload after damage is entirely directed to the intact columns that are closer to the damage area, 
and collapse propagates by nearest neighbor interactions. On the other hand, high stiffness of the beams 
and large plastic capacity of the columns induce democratic redistribution of overload between the columns. 
Consequently, the columns crush simultaneously triggering global pancake collapse. The collapse dynamics 

















































































































































































A= 0.75 
A= 1.33 


• brittle q^g 
° plastic qeq 


Figure 6: Equivalent collapse loads and residual strength fraction i?i for frames that undergo bending 
progressive collapse. 



He (= n+1) 


Figure 7: Collapse loads in Fig. 6 divided by L, considering the 2D frames as part of regular 3D structures. 


recorded in our simulations resembles global pancake. Note that in the studied framed structures, all the 
columns have identical compressive strength without disorder. Therefore, once the first two columns crush, 
pancake collapse can not be arrested. Nevertheless, at some fc > 0.35N/mm^ our frames undergo partial 
collapse because the progressive failure of the columns can be arrested by the initiation of bending collapse. 

Fig.m where superscript C indicates pancake collapse, shows that the collapse loads increase with the 
structural slenderness A = H/L, because the columns have tributary area related to and compressive 
strength proportional to H^. Furthermore, hierarchical structures with small n appear to be stronger than 
homogeneous ones both in terms of and of Qc- Finally, the residual strength fraction Ri toward pancake 
collapse is remarkably higher than that toward bending collapse (cf. Fig. |^, and is neither influenced by the 
hierarchical level 1/n, nor by A. In fact, Ri toward global pancake mode is related to the fraction of columns 
that are initially removed at one story. In our simulations, we always remove one third of the columns at one 
story, and obtain the constant value i?i « 0.6, slightly smaller than a theoretical 2/3 because of dynamics. 


Conclusions 

We showed how the dynamic strength after damage of 2D frames depends on the activated collapse 
mechanism. We can now drive a series of conclusions regarding the effect of damage extent, structural 


9 














t[s] 0.5 


1.5 


2.3 



> 


m. 


0.6 


0.9 


1.5 


4/3 


Figure 8: Snapshots of pancake collapse after damage for frames with n=5, A = 0.75, and n=ll, A = 1.33. 
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Figure 9: Equivalent loads and residual strength fraction for the studied structures undergoing pancake 
collapse. 


slenderness, and topological and mechanical hierarchy. Bending collapse provokes a local intensification 
of bending moments at the connections between the transfer beams above the damage area and the first 
intact column. Consequently, Qc and the residual strength fraction Ri decrease with the number of removed 
columns. In analogy with fracture mechanics, structures that are prone to bending collapse correspond to 
notch sensitive materials, and the number of removed columns corresponds to the crack width Chiaia ar 
Masoero 2008 . If global pancake collapse is triggered, Qc and Ri decrease with the fraction of removed 
columns, which is analogous to plastic failure of materials that are not notch sensitive. Consistently, Ri 
corresponding to global pancake collapse is remarkably larger than that corresponding to bending collapse. 

The structural slenderness X = H/L affects in general the collapse loads for both bending and pancake 
collapse modes. The effect of A depends on the scaling of cross section and reinforcement of the structural 


elements, with the beam length L and with the column height H (see the analytical results in Masoero 2010 


Masoero et al. 2013 , regarding the simulations in this paper). Nevertheless Ri turns out to be independent 


from A, because i?i is the ratio between two collapse loads with same scaling respect to L and H. 

Considering structural topology, in case of bending collapse hierarchical structures are more robust toward 
initial damage with fixed spatial extent (e.g. explosion, impact), and as robust as homogeneous structures 
toward single column removal (e.g. design error). The reason is that Ri toward bending collapse decreases 
with the number of removed columns at one story. This confirms the analogy with fracture mechanics, where 
notch sensitive hierarchical materials are tougher than homogeneous ones Lakes 1993 . On the other hand, 
considering global pancake collapse, structural hierarchy does not influence Ri toward initial damage with 
fixed spatial extent, while hierarchical structures are more sensitive than homogeneous ones to single column 
removals. This is due to i?i toward global pancake collapse decreasing with the fraction of removed columns. 
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FigJ shows that damaged frames undergoing global pancake collapse can carry the Ri « 60% of the static 
ultimate load of the intact structure Since well designed structures can carry a remarkably greater 
than the environmental load expected when an accidental event occurs, i?i related to global pancake collapse 
can ensure structural robustness for most of the practical cases Masoero 2010 . On the other hand, Ri 


related to bending collapse is generally much smaller, making structures vulnerable to accidental damage. 

In this work we considered idealized structures, with simplihed geometry and mechanical behavior of the 
elements. Reducing local complexity enables a better interpretation of the coral system response to damage. 
This study provides a basis of knowledge preceding the incorporation of more details and degrees of freedom, 
to investigate further aspects of progressive collapse. Shear failures can cause brittle ruptures and reduce 
the collapse resistance of large structural elements. Different locations of the initial damage may activate 
different collapse mechanisms. For example, damaging the upper stories would cause debris impacts, while 
removing external columns reduces qc without producing significant lateral toppling Calvi 2010 . The DEM 


algorithm was already applied to 3D structures in (Masoero, Wittel et ah, 2010), showing that the bending 
and pancake collapse mechanisms persist also in 3D. On the other hand, in 3D structures the horizontal 
floor slabs improve the horizontal redistribution of loads and the catenary action, increasing the strength 
toward bending collapse and impacting debris (see the Appendix and, e.g. , Vlassis et al. 2008 ). It is worth 


noting that horizontal ties and diaphragms increase the strength both after and before damage, causing a 
compensation that limits the effect on Ri. Finally, future works can incorporate a detailed description of 
structural connections, which are crucial for energy dissipation, catenary effect, and compartmentalization. 

Coming back to the central theme of structural hierarchy, our results already suggest that hierarchical 
structures are more robust toward accidental damage. An optimal solution would be to design: 1) a primary 
frame made of few large elements, with columns weaker than the beams, and 2) a secondary structure, 
made of many smaller elements, which defines the living space and follows traditional design rules. The 
primary frame would provide topological hierarchy, maximizing i?i toward bending collapse and enabling 
new possible compartmentalization strategies. The strong beams and weak columns of the primary frame 
would favor pancake collapse over bending collapse, and improve the vertical compartmentalization of high- 
rise buildings against falling debris. On the other hand, in real structures, the beams generally fail before 
the columns, and imposing the opposite is expensive. Nevertheless, designing a strong-beam weak-column 
behavior only for the primary frame can signihcantly limit the extra cost. Hierarchical structures can be 
a novel and somehow counterintuitive feature of robustness-oriented capacity design. Planning structural 
hierarchy requires understanding the complex system response to local damage, and should drive the design 
process since the very beginning. By contrast, traditional design is focused on local resistance against 
ordinary actions, and considers robustness toward accidents only at the end. This generally leads to non- 
hierarchical structures with strong columns and poorly understood system behavior. In addition, anti-seismic 
capacity design requires plastic failure of the beams to precede columns rupture (see e.g. [Byfield 2004 ). 
Overcoming these contradictions is a challenge toward optimizing structures against exceptional events. 
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Appendix: experimental benchmark 


In this appendix, we compare the numerical predictions of our DEM model with the experimental observa- 
2008 . We also briefly discuss some effects of catenary actions on collapse resistance. 


tions in Yi et al. 


Experiments The experimental setup in Yi et ah] 2008| consists of a plane frame made of reinforced 
concrete (see Fig. [Toja)). Columns are square in section (200x200mm), beams are rectangular (200mm tall, 
100mm wide). Everywhere, the longitudinal reinforcement is symmetrically distributed within the cross 


section (4(()12 steel bars). The strength and ultimate strain of concrete and steel are specified in |Yi et ah' 


2008| , while the elastic moduli are not. The mid column at the first floor is replaced by jacks that provides 
an upward vertical force N. In the middle of the top floor, a servo-hydraulic actuator applies a constant 
downward vertical force F=109kN, to represent the self weight of upper stories. Initially N = F =109kN, 
and then it is progressively reduced to reproduce quasi-static column loss, until a bending mechanism triggers 
collapse (see Fig.[Io|)b)). During the experiments, the increasing values of the midspan inflection A is plotted 
against A, to get the force-displacement reaction curve. The integral of the curve represents the energy 
dissipation capacity, which relates to the dynamic strength of the structure with respect to the activated 
collapse mechanism. 


Model description Our target is to capture the experimental reaction curve N — A through DEM simu¬ 


lations. The parametrization of our model, based on the geometry and mechanical data in Yi et al. 2008 


is straightforward. Therefore, we focus on the discrepancies between model and experimental inputs, and 
a few necessary additional assumptions. Regarding the overall geometry, we consider all the columns to be 
equally tall (1,100mm), while in the experiments the columns at the first floor were taller (1,567mm). This 
discrepancy should not have a significant effect on the collapse mechanism and the strength. The mechanical 
behavior of the real steel bars was strain hardening, with yielding at 416MPa, and rupture at 526MPa. In 
our model, we consider two limit cases of elastic-perfectly plastic behavior of the steel bars: weak steel “WS” 
with yielding threshold set at 416MPa, and strong steel “SS” yielding at 526MPa. Yi et al. 2008 provide 


two measures of the ultimate tensile strain S of the steel bars. We employ dio = 23%, which was measured 
on a longer bar segment, because in our simulation the strain develops within relatively long Euler-Bernoulli 
Beam Elements, EBEs. We assume Young moduli Eg = 200GPa for the steel, and = 30GPa for the 
concrete. In order to better understand the development of catenary actions, we consider two limit cases of 
cross section behavior under tension: fully reacting sections “FRS”, where the concrete always contributes 
to the tensile stiffness, and partially reacting sections “PRS” , where the concrete cracks and only the steel 
provides axial stiffness as soon as the cross section goes in tension. Furthermore, in order to focus exclusively 
on ruptures due to tensile strain in the steel, we allow for an infinite rotation capacity of the cross sections. 


Simulations and results We subject our model frames to gravity, but remain in the quasi-static regime 
by adding a high viscous damping force proportional to the velocity of each Spherical Discrete Element. 
We repeat numerous simulations with fixed F = I09kN and N, ranging from N = I09kN to values that 
are small enough to cause the quasi-static rupture of at least one EBE. We track the midspan deflection 
A (A)for comparison in Fig. [^c). 

Discussion In the experimental results, as N decreases, the system crosses several stages: (I) linear elastic 
A ^ 5mm, (II) elasto-plastic A < 22mm, (III) plastic hinges A < 140mm, (IV) catenary action A < 450mm, 
and (V) collapse. The transition from elastic to elasto-plastic is not evident from the curve, as well as that 
from plastic hinge to catenary action. By contrast, plastic hinges formation is clearly marked by a sudden 
change of slope at A « 40mm. Our simulations do not capture the initial elasto-plastic stage because we do 
not model the non-linear elasto-plastic behavior of concrete. This leads to an overestimation of the stiffness 
dN/dA before the formation of the plastic hinges. Nevertheless, the additional strain energy produced by 
this approximation is negligible when compared to the energy dissipated in the subsequent stages, ie. the 
overestimation of the initial stiffness is irrelevant for the actual dynamic collapse. Assuming weak steel WS, 
yielding at 4I6MPa, provides a good agreement with the experiment in terms of transition point to the plastic 
hinges stage. Gonsidering fractured concrete under tension yields the PRS-WS curve, which underestimates 
the structural strength at large A. The reason for this divergence can be that the steel hardens under strain. 
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Figure 10: (a) Geometry of the model frame and (b) generic bending collapse mechanism, (c) Experimental 
reaction curve, and predictions from our model. All the predicted curves are plotted until a first EBE 
fails in tension, with exception for the FRS-SS curve, truncated earlier due to evident divergence from the 
experimental result. If fully reactive cross sections FRS are assumed, collapse occurs when N << 0, i.e. N 
contributes to pulling down the structure instead of contrasting F. 


with reaction stress increasing from from 416MPa (WS) to 526MPa (SS). This interpretation is supported 
by the fact that the PRS-WS and PRS-SS curves envelop the experimental one. In particular, the PRS-SS 
curve reproduces well the last part of the experimental curve, as well as the collapse point. 

Conclusions Despite strong simplifying assumptions in the formulation, our DEM model provides rea¬ 
sonably good quantitative predictions of the experimental results. For the simulations in the body of this 
paper, we always considered fully reactive cross sections with concrete that does not crack under tension. 
The FRS-WS curve in Fig. [To)[c) shows that this assumption leads to an overestimation of the static collapse 
strength against column removal (« -1-70%). Let us conjecture that FRS induce the same strength increase 
of -|-70% in the structure without column removal, i.e. in q^. From a heuristic application of energy conser¬ 
vation, one can estimate the dynamic collapse load after sudden column removal by considering the mean 
N in the catenary stage: N = 20kN from the experiment, and N = — 20kN from the simulation with FRS- 
WS. Consequently, the increase in post-damage dynamic collapse strength qc due to FRS is approximately 
(1 — 130/90) • 100 = -1-44%. In conclusion, assuming fully reactive cross sections causes an underestimation 
of the residual strength fraction Ri = qc/qi, which our example quantifies as (44/70 — 1) « —37%. However 
this assumption does not affect the main statement of this work on hierarchical structures. 
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